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1 . INTRODUCTION 


There have been efforts to improve the color film product that is currently 
being used in the Large Area Crop Inventory Experiment (LACIE). Most of 
those efforts (refs. 1, 2) focused on the gain and bias factors used prior - 
to the generation of imagery. In this study, a method is developed to perform 
a linear transformation from the Land Sate^lille (Landsat) data space into the 
color gun space of the film converter so that the data will occupy a maximum 
volume in a cube defined in the color gun space. Limitations on the volume 
expansion are expressed in terms of placing gun saturation a specified number 
of standard deviations from the mean and limiting the expansion ratio so that 
noise-level variations in the data are not quite visible to the eye. Landsat 
data bright,hess (using the first row in the Kauth transformation, ref. 3) is 
aligned with the gray axis of the blue (B), green (G), and red (R) gun space 
of the production film converter (PFC). The Kauth greenness and yellowness 
are presented orthogonally to each other and to the brightness in the color 
gun space. 

A description of the transformation involved in this new color film generation 
method is presented in section 2. Programming mechanics and examples are 
given in section 3. A subjective evaluation of this new product, along with 
some remarks, is presented in section 4. 

2. DESCRIPTION OF THE TRANSFORMATION 

Let Xj^, Xg, and x^ be Kauth brightness, greenness, and yellowness, respec- 
tively. Obtain a, $, and y as follows: 
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where (x. , a. ), {x , o ), and (x , a ) are the mean and standard deviation of 
u M y y y j 

the Kauth brightness, greenness, and yellowness, respectively, and n is the 
number of standard deviations to be used for inclusion of the data in the 
gun cube. Equation (1) amounts to standardizing the brightness, greenness, 
and yellowness, then applying a gain to the brightness so that (xj^ ± noj^) 

255i/3’ 

falls within ± — 2 — » while applying a gain to the greenness and yellowness 

so that (x ± na ) and (x ± na ) falls within ±-^| — . Note that the gray 
axis of the PFC color gun cube is of length 255/3. 

Let the and ^ color guns be three Cartesian coordinates in Euclidean 

space, where the bar under G, B, and R is used to indicate a vector. In 
equation (1), £ is aligned with the gray axis of the color gun cube, i.e., 

It points in the direction from (0, 0, 0) to (255, 255, 255) in the space 
defined above. Note that (1//3) a is the component along each of the 6, 
and ^ color guns. (Notice that since ^ can take on negative values, the 
components of along each color gun can be negative. This will be compen- 
sated for later by adding a bias term.) Looking into the gray axis of the 

color gun cube as shown in figure 1, £ is defined as lying in the {^, a^) 
plane with ^ • a = 0 and • R > 0, where “•*' indicates the dot product. 

Def 1 ne 11 as : 

X a (2) 

where "x" indicates the cross product. Angle 6 in figure 1 is the remaining 
degree of freedom with 

?_’= (y cos 0 - b sin (3) 

jl = (B cos 0 + Y sin 0 ) 11 ^. (4) 

where ^ and jr. are unit vectors along the directions of ^ and n_, respectively. 
Define 1 , j_, and j< as unit vectors along the and R vectors, respectively. 

Resolving and ji into components along the G, and R directions 

(fig. 2) yields: 

^ = a cos ^(i_ + i + Ji) 

£ = I sin <{• li + ^ 
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(5) 

( 6 ) 




Figure 1.- Looking into the gray 


i • a = 0 

Y * g. = 0 

Y, 3_, and are all 
perpenJicular to a. 

C • Y “ cos 0 = n • 3 


axis of the color gun cube. 
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cos (}>(i + i)cos 45° (7) 

rt^ = n cos 45 ° (t_ - (8) 

where cos 4« = and sin <|) = /2/3. Each of the color gun axes is the sum 

of its components; hence, 

R = a cos ^ ? sin ^ (9) 

G = a cos - C cos <(> cos 45° + n cos 45° (10) 

B = a cos - C cos <f> cos 45° - n cos 45° (11) 


By substituting equations (3) and (4) for 5 and n m equations (10) and (11) 
and then adding a bias to center the color gun distributions in the cube 
one obtains* 


n 1 — ^ 

Vf 


- sin 6 e + cos 9 Y + 


255 


G‘ 

B' 


_ a ^ /sin e ^ cos e\ /sin e cos eV. , 255 

■ ^ ~r ~r ~ 

sin e\ , 255 

jz r 2 


ct , /sin e cos e 
^ 


\ /cos 0 

r U 


3. PROGRAMMING MECHANICS AND EXAMPLES 


(12) 

(13) 

(14) 


In equations (12), (13), and (14), it is desirable to compute the R', B', and 
6' color gun values using a method that will quantize the color gun levels 
only once; i.e. , just prior to storing those levels on a magnetic tape. In 
fact, once ct, 6, and y are computed, it is not possible to store them on a 
universally formatted magnetic tape because a, 3, and y range from negative 
to positive values. The method used to calculate B‘, G', and R' is described 
as follows. 
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Let C = (B'. 6‘. = (255/2, 255/2, 255/2}''^, v = (a, B, y)'^ , where T 

IS the transpose. Define the 3 by 3 matrix M as: 


M = 


(1//3) (sin 6/t/S - cos q/}/2) -(cos e//6 + sin 9//2) 

(1//5) (sin 6/ /S’ + cos Q/t/2) (sin Q//Z - cos Q/^/6) 

.(1//3) -v^ sin e cos e 


' <™ia) 

Equations (14), (13), and (12) , .respectively , can be written as: 


(15) 


£ = My^ + m;j 


(16) 


Let ^ be a four-dimensional feature vector from Landsat data with mean m and 

covariance matrix s . Let K be a (3 x 4) matrix formed from the first three 

^ 2 2 2 
.rows of the Kauth transformation matrix. Let cr. , u , and c be the first, 

b g’ y 

T 

second, and third diagonal elements of the matrix Kz K , respectively. 

Define the matrix A as: 


■'255/3/2na, 


A = 


255v^/2n0 


0 


255v^/2na 


The vector y^ in equation (16) can be written as: 

V = AK^ - AKm 

Substituting equation (18) into equation (16) yields: 

C = AX + b 


where 


(17) 


(18) 

(19) 


A = MAK 
£ = nij - MAKm 


(20) 

(21) 


Equation (19) was applied to nine LAC IE segments with each segment having 
four acquisitions in the 1976 crop year. Some examples of the new product 
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as well as the current LACIE product are shown in figure 3. The new product 
shown in figure 3 was generated for n = 3 and 0 = -n/2. The value of e was 
chosen to be “ir/2 in order to align the greenness mostly with the red color 
gun of the PFC. In comparison with the current LACIE product, the new 
product has more contrast and colors and appears to be brighter* The field 
boundaries in the new product are more pronounced than in the current LACIE 
product. 


4. REMARKS 

An image analyst (private communication with W. T. Hocutt of Lockheed Elec- 
tronics Company, Inc., Systems and Services Division) performed a subjective 
evaluation of the new product and reported the following. 

a. Most acquisitions showed speckled appearance. This makes field identifi- 
cation difficult. 

b. The purple signature in the new product is difficult to define. 

c. in segment 1645, fields that are similar on the current LACIE product are 
different on the new product, and fields that are different on the current 
LACIE product are similar on the new product. The speckled appearance 

in the new product is primarily due to excessive stretching in the yellow- 

ness direction. The yellowness variance of all acquisitions was 

examined and found to be between 1.31 and 14.66. A 14.66 yellowness 

variance corresponded to about a 2-percent cloud covered acquisition. 

2 

An acceptable level of speckle was found to correspond to = 5.94. 

Hence, one way of removing most of the speckle is by multiplying y by a 
factor F, where 

Z 

F = i (22) 

Max(a , 6.0) 

*J 

This product produces maximal chromatic expansion of data, with certain con- 
straints. The product is not a "true color" image, two points in data space 
that are linearly related do not, in general,' share such chromatic attributes 
as hue and saturation when mapped into color space under this transformation. 
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(a) Current LACIE color film product of segment 1618 in Grand Forks, 
North Dakota. Acquisition date: August 22, 1976. 



(b) New color film product of segment 1618. 


Figure 3.- Some examples of the LACIE current and new color film products. 
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(c) Current LACIE color film product of segment 15^1. 
Acquisition date: July 23, 1976. 



(d) Mew color film product of segment 1541. 


Figure 3.— Concluded. 
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